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(54) Magnetic field in homogeneity measurement method and apparatus, Phase correction 
method and apparatus 

(57) In order to perform imaging with water and fat 
separated in an efficient manner, a tomographic image 
is obtained having a phase difference of n/2 between 
water and fat using magnetic resonance, the phase of 
the tomographic image is multiplied by n, the phase 
after being corrected for wraparound is multiplied by 
1/n, and a phase image representing magnetic field 
inhomogeneity is generated. The tomographic image is 
corrected for phase using the phase image, and is sep- 
arated into a water image and a fat image. 
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Description 

[0001] The present invention relates to a magnetic 
field inhomogene'rty measurement method and appara- 
tus, a phase correction method and apparatus and a 5 
magnetic resonance imaging apparatus, and more par- 
ticularly to a method and apparatus for measuring mag- 
netic field inhomogeneity in a space to be imaged based 
on the phase of pixel data of an image obtained by mag- 
netic resonance imaging, a method and apparatus for 10 
correcting the phase of pixel data based on the meas- 
ured magnetic field inhomogeneity, and a magnetic res- 
onance imaging apparatus for obtaining images with 
water and fat separated based on the phase-corrected 
pixel data. is 
[0002] A magnetic resonance imaging apparatus is 
configured to generate a static magnetic field in a space 
containing a subject to be imaged, generate gradient 
and high frequency magnetic fields in the static mag- 
netic field space, and produce (i.e., reconstruct) an 20 
image based on magnetic resonance signals generated 
from the magnetic field space. Since the magnetic reso- 
nance signals from fat have a different frequency than 
the magnetic resonance signals from water due to 
chemical shift, water and fat can be separately imaged 25 
utilizing a phase difference based on the difference in 
frequency. 

[0003] The phase of the magnetic resonance sig- 
nals is affected by inhomogeneity of the static magnetic 
field strength, and imaging by a so-called Dixon method 30 
is performed to separately image water and fat without 
being affected by the magnetic field inhomogeneity. 
This technique involves obtaining an image with the 
phase of magnetic resonance signals from water (which 
will be simply referred to as a water phase he rein be low) 35 
in phase with that of magnetic resonance signals from 
fat (which will be simply referred to as a fat phase here- 
inbelow) and an image with the water phase opposite to 
the fat phase, and generate a water image through the 
sum of these images and a fat image through the differ- 40 
ence between these images. 

[0004] However, the Dixon method requires normal 
imaging operations for two images to obtain one image, 
resulting in low efficiency. Moreover, since the phase 
change rate for fat decreases relative to that for water as 45 
the static magnetic field strength lowers, a longer TE 
(echo time) is required to obtain signals with water and 
fat in phase by using, for example, a gradient echo, and 
accordingly signal decay increases, making it impracti- 
cal to employ the Dixon method in a magnetic field hav- 50 
ing a strength as low as 0.2 T, for example. 
[0005] The present invention therefore seeks to 
provide a magnetic field inhomogeneity measurement 
method and apparatus for performing imaging with 
water and fat separated in an efficient manner, a phase 55 
correction method and apparatus based on the mag- 
netic field inhomogeneity, and a magnetic resonance 
imaging apparatus for performing imaging with water 



and fat separated in an efficient manner. 
[0006] In accordance with a first aspect of the 
invention, there is provided a magnetic field inhomoge- 
neity measurement method comprising the steps of: 
obtaining a tomographic image of a subject within a 
static magnetic field space using magnetic resonance, 
the tomographic image having a phase difference of 
2n/n (n > 2) between pixel data of water and fat; multi- 
plying the phase of pixel data of the tomographic image 
by n; correcting wraparound caused by the phase multi- 
plied by n exceeding a range between ±71; and multiply- 
ing the wraparound-corrected phase by 1/n to obtain a 
phase representing magnetic field inhomogeneity in the 
static magnetic field space. 

[0007] In accordance with a second aspect of the 
invention, there is provided a magnetic field inhomoge- 
neity measurement apparatus comprising: imaging 
means for obtaining a tomographic image of a subject 
within a static magnetic field space using magnetic res- 
onance, the tomographic image having a phase differ- 
ence of 2n/n (n > 2) between pixel data of water and fat; 
phase-multiplier means for multiplying the phase of 
pixel data of the tomographic image by n; wraparound- 
correction means for correcting wraparound caused by 
the phase multiplied by n exceeding a range between 
±ji; and magnetic field in homogeneity-measurement 
means for multiplying the wraparound-corrected phase 
by 1/n to obtain a phase representing magnetic field 
inhomogeneity in the static magnetic field space. 
[0008] In accordance with a third aspect of the 
invention, there is provided a phase correction method 
comprising the steps of: obtaining a tomographic image 
of a subject within a static magnetic field space using 
magnetic resonance, the tomographic image having a 
phase difference of 2n/n (n > 2) between pixel data of 
water and fat; multiplying the phase of pixel data of the 
tomographic image by n; correcting wraparound caused 
by the phase multiplied by n exceeding a range between 
±tc; multiplying the wraparound-corrected phase by 1/n 
to obtain a phase representing magnetic field inhomo- 
geneity in the static magnetic field space; and correcting 
the phase of pixel data of the tomographic image using 
the phase obtained. 

[0009] In accordance with a fourth aspect of the 
invention, there is provided a phase correction appara- 
tus comprising: imaging means for obtaining a tomo- 
graphic image of a subject within a static magnetic field 
space using magnetic resonance, the tomographic 
image having a phase difference of 2rc/n (n > 2) between 
pixel data of water and fat; phase-multiplier means for 
multiplying the phase of pixel data of the tomographic 
image by n; wraparound-correction means for correct- 
ing wraparound caused by the phase multiplied by n 
exceeding a range between ±n; magnetic field inhomo- 
geneity-measurement means for multiplying the wrapa- 
round-corrected phase by 1/n to obtain a phase 
representing magnetic field inhomogeneity in the static 
magnetic field space; and phase-correction means for 
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correcting the phase of pixel data of the tomographic 
image using the phase obtained. 

[0010] In accordance with a fifth aspect of the 
invention, there is provided a magnetic resonance imag- 
ing apparatus comprising: imaging means for obtaining 
a tomographic image of a subject within a static mag- 
netic field space using magnetic resonance, the tomo- 
graphic image having a phase difference of 2n/n (n > 2) 
between pixel data of water and fat; phase-multiplier 
means for multiplying the phase of pixel data of the tom- 
ographic image by n; wraparound-correction means for 
correcting wraparound caused by the phase multiplied 
by n exceeding a range between ±n; magnetic field 
inhomogeneity-measurement means for multiplying the 
wraparound-corrected phase by 1/n to obtain a phase 
representing magnetic field inhomogeneity in the static 
magnetic field space; phase-correction means for cor- 
recting the phase of pixel data of the tomographic image 
using the phase obtained; and image-producing means 
for separately producing a water image and a fat image 
using the phase difference in the phase-corrected pixel 
data. 

[0011] In accordance with a sixth aspect of the 
invention, there is provided a magnetic resonance imag- 
ing method comprising the steps of: obtaining a tomo- 
graphic image of a subject within a static magnetic field 
space using magnetic resonance, the tomographic 
image having a phase difference of 2x/n (n > 2) between 
pixel data of water and fat; multiplying the phase of pixel 
data of the tomographic image by n; correcting wrapa- 
round caused by the phase multiplied by n exceeding a 
range between ±k\ multiplying the wraparound-cor- 
rected phase by 1/n to obtain a phase representing 
magnetic field inhomogeneity in the static magnetic field 
space; correcting the phase of pixel data of the tomo- 
graphic image using the phase obtained; and sepa- 
rately producing a water image and a fat image using 
the phase difference in the phase-corrected pixel data. 
[0012] In accordance with a seventh aspect of the 
invention, there is provided a magnetic resonance imag- 
ing method comprising the steps of: obtaining a tomo- 
graphic image of a subject within a static magnetic field 
space using magnetic resonance, said tomographic 
image having a phase difference of 2n/n (n > 2) between 
pixel data of water and fat; multiplying the phase of pixel 
data of said tomographic image by n; correcting wrapa- 
round caused by said phase multiplied by n exceeding a 
range between ±rc; multiplying said wraparound-cor- 
rected phase by 1/n to obtain a phase representing 
magnetic field inhomogeneity in said static magnetic 
field space; correcting the phase of pixel data of said 
tomographic image using said phase obtained; and 
separately producing a water image and a fat image 
using the phase difference in said phase-corrected pixel 
data. 

[0013] In this latter aspect, the method may com- 
prise the steps of: modifying the phase difference in 
said phase-corrected pixel data based on a histogram of 



the phase of said pixel data; and separately producing a 
water image and a fat image using said modified phase 
difference. 

[0014] In all aspects, n may be equal to 4. 

5 [0015] Thus, with the present invention, water and 
fat are made in phase by multiplying the phase of pixel 
data by n, and wraparound that exceeds a range 
between ±n is corrected and then is multiplied by 1/n to 
obtain a phase distribution corresponding to magnetic 

10 field inhomogeneity. The phase distribution is used to 
correct the phase of pixel data, thereby eliminating the 
effect of the magnetic field inhomogeneity. Using the 
pixel data from which the effect of the magnetic field 
inhomogeneity has been eliminated, separate water 

15 and fat images are provided using the phase difference. 
[0016] Thus, the present invention can provide a 
magnetic field inhomogeneity measurement method 
and apparatus for performing imaging with water and fat 
separated in an efficient manner, a phase correction 

20 method and apparatus based on the magnetic field 
inhomogeneity, and a magnetic resonance imaging 
apparatus for performing imaging with water and fat 
separated in an efficient manner. 
[0017] The invention will now be described in 

25 greater detail, by way of example, with reference to the 
drawings, in which: 

Rgure 1 is a block diagram of an apparatus -in 
accordance with an embodiment of the present 
30 invention. 

Rgure 2 is a block diagram of an apparatus in 
accordance with another embodiment of the 
present invention. 

35 

Rgure 3 illustrates an exemplary pulse sequence 
executed by the apparatus shown in Fig. 1 or 2. 

Rgure 4 illustrates another exemplary pulse 
40 sequence executed by the apparatus shown in Rg. 
1 or 2. 

Rgure 5 is a block diagram of a data processing 
section in the apparatus shown in Rg. 1 or 2. 

45 

Rgures 6 — 9 illustrate the function of the data 
processing section shown in Rg. 5. 

Rgure 10 is a graphic representation illustrating a 
so histogram of the phase of a complex image. 

Rgure 1 1 is a vector diagram of water and fat sig- 
nals. 

55 [0018] Figure 1 shows a block diagram of a mag- 
netic resonance imaging apparatus, which is an embod- 
iment of the present invention. The configuration of the 
apparatus represents an embodiment of the apparatus 
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in accordance with the present invention, and the oper- 
ation of the apparatus represents an embodiment of the 
method in accordance with the present invention. 
[0019] As shown in Fig. 1, the present apparatus 
has a magnet system 100. The magnet system 1 00 has 
a main magnetic field coil section 102, a gradient cot! 
section 106 and an RF (radio frequency) coil section 
108. These coil sections have a generally cylindrical 
outer shape and are concentrically disposed. A subject 
to be imaged 300 is rested on a cradle 500 and carried 
into and out of the internal space of the magnet system 
100 by carrier means (not shown). The subject to be 
imaged 300 is an embodiment of the subject of the 
present invention. 

[0020] The main magnetic field coil section 102 
generates a static magnetic field in the internal space of 
the magnet system 1 00. The direction of the static mag- 
netic field is generally in parallel with the direction of the 
body axis of the subject to be imaged 300, i.e., a so- 
called horizontal magnetic field is generated. The main 
magnetic field coil section 1 02 comprises a supercon- 
ductive coil, for example. It will be easily recognized that 
the main magnetic field coil section 1 02 is not limited to 
comprising a superconductive coil, but may comprise a 
normal conductive coil. 

[0021] The gradient coil section 106 generates gra- 
dient magnetic fields for imparting gradients to the static 
magnet field strength. The gradient magnetic fields to 
be generated include three types of magnetic fields: a 
slice gradient magnetic field, a readout gradient mag- 
netic field and a phase encoding gradient magnetic 
field, and the gradient coil section 106 has three gradi- 
ent coils (not shown) corresponding to these three gra- 
dient magnetic fields. 

[0022] The RF coil section 108 generates a high 
frequency magnetic field in the static magnetic field 
space to excite spins within the subject to be imaged 
300. The generation of the high frequency magnetic 
field will be referred to as transmission of an RF excita- 
tion signal hereinbelow. The RF coil section 108 also 
receives an electromagnetic wave generated by the 
excited spins, i.e., a magnetic resonance signal. The RF 
coil section 1 08 has a transmit coil and a receive coil 
(not shown). The transmit and receive coils may be 
either the same coil or separate dedicated coils. 
[0023] The gradient coil section 106 is connected 
with a gradient driving section 1 30 for supplying driving 
signals to the gradient coil section 106 to generate the 
gradient magnetic fields. The gradient driving section 
130 has three driving circuits (not shown) correspond- 
ing to the three gradient coils in the gradient coil section 
106. 

[0024] The RF coil section 1 08 is connected with an 
RF driving section 140 for supplying driving signals to 
the RF coil section 1 08 to transmit the RF excitation sig- 
nals, thereby exciting the spins within the subject to be 
imaged 300. The RF coil section 108 is also connected 
with a data acquisition section 150 for gathering signals 



received by the RF coil section 108 and collecting the 
signals as digital data. 

[0025] The gradient driving section 130, RF driving 
section 140 and data acquisition section 150 are con- 
5 nected with a control section 160 for controlling these 
sections 130 — 150. 

[0026] The output of the data acquisition section 
150 is connected to a data processing section 170. The 
data processing section 1 70 stores data gathered from 

10 the data acquisition section 150 in a memory (not 
shown). Thus, a data space is formed in the memory, 
which constitutes a two-dimensional Fourier space. The 
data processing section 170 performs an inverse two- 
dimensional Fourier transformation on the data in the 

15 two-dimensional Fourier space to reconstruct an image 
of the subject to be imaged 300. 

[0027] The data processing section 170 is con- 
nected to the control section 160. The data processing 
section 170 is superior to the control section 160 and 
20 governs the section 160. The data processing section 
170 is connected with a display section 180 that dis- 
plays the reconstructed image and several information 
output from the data processing section 170, and an 
operating section 1 90 that is operated by a human oper- 
as ator and inputs several commands, information and so 
forth to the data processing section 170. 
[0028] Figure 2 shows a block diagram of another 
magnetic resonance imaging apparatus, which is an 
embodiment of the present invention. The configuration 
30 of the apparatus represents an embodiment of the 
apparatus in accordance with the present invention, and 
the operation of the apparatus represents an embodi- 
ment of the method in accordance with the present 
invention. 

35 [0029] The apparatus shown in Fig. 2 has a differ- 
ent magnet system 100' than the apparatus shown in 
Fig. 1 . Components except the magnet system 100' are 
configured in a similar manner to those in the apparatus 
shown in Fig. 1, and similar parts are designated by 

40 similar reference numerals and explanation thereof will 
be omitted. 

[0030] The magnet system 100' has a main mag- 
netic field magnet section 102 1 , a gradient coil section 
106' and an RF coil section 108'. The main magnetic 

45 field magnet section 102' and the coil sections 106' and 
1 08' each comprise a pair of respective members facing 
each other with a space interposed. These sections 
have a generally disk-like outer shape and are disposed 
with a common center axis. The subject to be imaged 

so 300 is rested on the cradle 500 and carried into and out 
of the internal space of the magnet system 1 00' by car- 
rier means (not shown). 

[0031 ] The main magnetic field magnet section 1 02' 
generates a static magnetic field in the internal space of 
55 the magnet system 1 00'. The direction of the static mag- 
netic field is generally orthogonal to the direction of the 
body axis of the subject to be imaged 300, i.e., a so- 
called vertical magnetic field is generated. The main 
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magnetic field magnet section 1 02' comprises a perma- 
nent magnet, for example. It will be easily recognized 
that the main magnetic field magnet section 102* is not 
limited to comprising a permanent magnet, but may 
comprise a superconductive or normal conductive elec- 
tromagnet, etc. 

[0032] The gradient coil section 1 06' generates gra- 
dient magnetic fields for imparting gradients to the static 
magnet field strength. The gradient magnetic fields to 
be generated include three types of magnetic fields: a 
slice gradient magnetic field, a readout gradient mag- 
netic field and a phase encoding gradient magnetic 
field, and the gradient coil section 106' has three gradi- 
ent coils (not shown) corresponding to these three gra- 
dient magnetic fields. 

[0033] The RF coil section 108' transmits an RF 
excitation signal to the static magnetic field space to 
excite spins within the subject to be imaged 300. The 
RF coil section 108 also receives a magnetic resonance 
signal generated by the excited spins. The RE coil sec- 
tion 108' has a transmit coil and a receive coil (not 
shown). The transmit and receive coils may be either 
the same coil or separate dedicated coils. 
[0034] Figure 3 shows an exemplary pulse 
sequence for use in magnetic resonance imaging. The 
illustrated pulse sequence is one for a spin echo (SE) 
method. 

[0035] Specifically, (1) is a sequence of 90° and 
180° pulses for RF excitation of the SE method, and (2), 
(3), (4) and (5) are sequences of a slice gradient Gs, a 
readout gradient Gr, a phase encoding gradient Gp and 
a spin echo MR, respectively, of the SE method. It 
should be noted that the 90° and 1 80° pulses are repre- 
sented by their respective central signals. The pulse 
sequence proceeds from the left to the right along a 
time axis t. 

[0036] As shown, the 90° pulse results in 90° exci- 
tation of the spins. At the same time, the slice gradient 
Gs is applied to achieve selective excitation for a certain 
slice. After a predetermined period from the 90° excita- 
tion, 180° excitation by the 180° pulse, or spin inversion, 
is performed. Again, the slice gradient Gs is applied at 
the same time to achieve selective inversion for the 
same slice. 

[0037] During the time between the 90° excitation 
and the spin inversion, the readout gradient Gr and 
phase encoding gradient Gp are applied. The readout 
gradient Gr dephases the spins, and the phase encod- 
ing gradient Gp phase-encodes the spins: 
[0038] After the spin inversion, the spins are 
rephased by the readout gradient Gr to generate a spin 
echo MR. The spin echo MR is collected by the data 
acquisition section 150 as view data. Such a pulse 
sequence is repeated 128 — 256 times in a cycle of TR 
(repetition time). The phase encoding gradient Gp is 
varied for each repetition to provide different phase 
encoding each time. Thus, view data for 128 — 256 
views are obtained. 



[0039] The spin echo MR is an RF signal having a 
symmetrical waveform with respect to the echo center. 
The central echo occurs after TE (echo time) from the 
90° excitation. By appropriately selecting the time TE, 

5 the phase difference between echoes of water and fat 
can be 2n/n (n > 2). For example, when n = 4, the value 
of TE for a phase difference of tU2 is of the order of 2x + 
8.6 ms or 2x - 8.6 ms under a static magnetic field 
strength of 0.2 T, wherein x is a time interval from 90° 

to excitation to 180° excitation. The spin echo obtained at 
such TE has a sufficient signal strength. 
[0040] Another example of the pulse sequence for 
magnetic resonance imaging is shown in Fig. 4. This 
pulse sequence is one for a gradient echo (GRE) 

15 method. 

[0041] Specifically, (1 ) is a sequence of an ct° pulse 
for RF excitation of the GRE method, and (2), (3), (4) 
and (5) are sequences of a slice gradient Gs, a readout 
gradient Gr, a phase encoding gradient Gp and a gradi- 
20 ent echo MR, respectively, of the GRE method. It should 
be noted that the ot° pulse is represented by its central 
signal. The pulse sequence proceeds from the left to the 
right along a time axis t. 

[0042] As shown, the cc° pulse results in cc° excita- 
25 tion of the spins, wherein a is not greater than 90. At the 
same time, the slice gradient Gs is applied to achieve 
selective excitation for a certain slice. 
[0043] After the ct° excitation, the spins are phase- 
encoded by the phase encoding gradient Gp. Next; the 
30 spins are first dephased and are subsequently 
rephased by the readout gradient Grto generate a gra- 
dient echo MR. The gradient echo MR is collected by 
the data acquisition section 150 as view data. Such a 
pulse sequence is repeated 128 — 256 times in a cycle 
35 of TR. The phase encoding gradient Gp is varied for 
each repetition to provide different phase encoding 
each time. Thus, view data for 128 — 256 views are 
obtained. 

[0044] The gradient echo MR is an RF signal hav- 
40 ing a symmetrical waveform with respect to the echo 
center. The central echo occurs after TE from the a° 
excitation. By appropriately selecting the time TE, the 
phase difference between echoes of water and fat can 
be 2ji/n (n > 2). For example, when n = 4, the value of 
45 TE for a phase difference of it/2 is of the order of 8.6 ms 
under a static magnetic field strength of 0.2 1 . The gra- 
dient echo obtained at such TE has a sufficient signal 
strength. 

[0045] The view data obtained by the pulse 
so sequence of Fig. 3 or 4 are collected into the memory in 
the data processing section 170. It will be easily recog- 
nized that the pulse sequence is not limited to the SE or 
GRE method, but may be of any other appropriate tech- 
nique such as a fast spin echo (FSE) method. 
55 [0046] The data processing section 170 performs 
an inverse two-dimensional Fourier transformation on 
the view data to reconstruct a tomographic image of the 
subject to be imaged 300. The reconstructed image is 
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stored in the memory. The configuration and function of 
the present apparatus in the preceding description rep- 
resent an embodiment of the imaging means of the 
present invention. 

[0047] The data processing section 170 produces 
an image rendering water and an image rendering fat 
from the reconstructed image. The image rendering 
water will be referred to as a water image, and the 
image rendering fat as a fat image hereinbelow. 
[0048] in producing the water and fat images, the 
data processing section 1 70 generates a strength distri- 
bution of the static magnetic field. The data processing 
section 170 is an embodiment of the magnetic field 
inhomogeneity measurement apparatus of the present 
invention. The configuration of the data processing sec- 
tion 170 represents an embodiment of the apparatus in 
accordance with the present invention, and the opera- 
tion thereof represents an embodiment of the method in 
accordance with the present invention. 
[0049] The data processing section 170 also per- 
forms phase correction for eliminating the effect of the 
magnetic field inhomogeneity based on the magnetic 
field inhomogeneity. The data processing section 170 is 
an embodiment of the phase correction apparatus of the 
present invention. The configuration of the data 
processing section 170 represents an embodiment of 
the apparatus in accordance with the present invention, 
and the operation thereof represents an embodiment of 
the method in accordance with the present invention. 
[0050] Figure 5 shows a block diagram of the data 
processing section 1 70 with respect to separate produc- 
tion of water and fat images. The function of the blocks 
in Fig. 5 is implemented by, for example, a computer 
program. 

[0051] As shown, the data processing section 170 
has a power image producing section 702 and a phase 
image producing section 704, which are supplied with a 
reconstructed image from a preceding image recon- 
struction section 700. Pixel data of the reconstructed 
image is represented by a complex number. That is, the 
pixel data has a real component and an imaginary com- 
ponent. The real component will be referred to as a real 
part and the imaginary component as an imaginary part 
hereinbelow. 

[0052] The power image producing section 702 
generates the power of complex data for each pixel to 
produce an image using the power as each pixel value, 
i.e., a power image. The phase image producing section 
704 generates the phase of complex data for each pixel, 
i.e., the arc tangent of the real part and the imaginary 
part, to produce an image using the phase as each pixel 
value, i.e., a phase image. 

[0053] Figure 6(a) shows a schematic of the phase 
image, in which a one-dimensional profile of the phase 
image is illustrated when the tomographic image com- 
prises a fat image and a water image surrounding the fat 
image. Figure 6(a) illustrates the case in which the 
phase difference between the water and fat images is 



7i/2, i.e., n is set to 4 for the phase difference of 27t/n. 
While the case of n = 4 will be described hereinbelow, 
the description also applies to other values of n. 
[0054] If the static magnetic field is homogeneous, 

5 the one-dimensional profile of the phase image (which 
will be simply referred to as a phase image) is to have a 
shape indicated by dot-dash line in Fig. 6(a) because 
the phase of water image is zero. However, if the static 
magnetic field has linearly slanting inhomogeneity, for 

10 example, the phase image will have a shape indicated 
by solid line. 

[0055] The phase image is input to a phase x n mul- 
tiplier section 706, which is an embodiment of the phase 
multiplier means of the present invention. The phase x n 

15 multiplier section 706 multiplies the pixel values of the 
phase image by n (= 4). Thus, a phase image as shown 
in Fig. 6(b) is obtained. As shown, the phase difference 
between water and fat becomes 2n by multiplication by 
4, causing water and fat to be in phase. In the phase 

20 image, so-called wraparound occurs in which a portion 
beyond the range ±k turns into the range ±n. In addition, 
discontinuity or abrupt change of the phase occurs at 
the boundary between water and fat. 
[0056] Such a phase image is input to a complex 

25 image producing section 708. The complex image pro- 
ducing section 708 is also supplied with the power 
image from the power image producing section 702. 
The complex image producing section 708 produces a 
complex image based on the phase image and the 

30 power image. 

[0057] The real part of the complex image can be 
obtained as the cosine of the power image data, and the 
imaginary part can be obtained as the sine of the power 
image data. It should be noted that the angle for use in 

35 the cosine and sine calculation is a phase angle. 

[0058] The complex image is input to a phase 
image producing section 712 through a low-pass filter 
section 710. The phase image producing section 712 
produces a phase image from the low-pass filtered 

40 complex image. In the phase image to be produced, the 
discontinuous or abrupt changing phase portion as 
shown in Fig. 7(a) is made continuous or alleviated by 
the low-pass filtering as exemplarily shown in Fig. 7(b). 
[0059] Such a phase image is input to a phase 

45 unwrapping section 71 4, which is an embodiment of the 
wraparound-correction means of the present invention. 
The phase unwrapping section 714 unwraps the wrap- 
ping phase shown in Fig. 8(a) into the phase as shown 
in Fig. 8(b). 

so [0060] The unwrapped phase image is input to a 
phase x 1/n multiplier section 716, which is an embodi- 
ment of the magnetic field inhomogeneity measurement 
means of the present invention. The phase x 1/n multi- 
plier section 716 multiplies the pixel values of the input 

55 phase image by 1/n (= 4). Thus, a phase image as 
shown in Fig. 8(c) is obtained. This phase image corre- 
sponds to a phase image of the subject to be imaged 
300 that consists only of water. Therefore, the phase 
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image represents the strength distribution of the static 
magnetic field, or static magnetic field in homogeneity. 
[0081] The phase image representing the static 
magnetic field inhomogeneity is input to a phase correc- 
tion section 71 8, which is an embodiment of the phase- 
correction means of the present invention. The phase 
correction section 718 is also supplied with the phase 
image from the phase image producing section 704, 
which is a phase image of water and fat containing a 
phase error due to the static magnetic field inhomoge- 
neity. The phase correction section 718 corrects the 
phase error in the phase image of water and fat using 
the phase image representing the static magnetic field 
inhomogeneity. 

[0062] Specifically, as shown in Fig. 9, the pixels in 
the phase image representing the static magnetic field 
inhomogeneity (b) is subtracted from the corresponding 
pixels in the phase image of water and fat containing the 
phase error (a) to generate a phase image not contain- 
ing the phase error due to the static magnetic field inho- 
mogeneity as shown in (c). 

[0063] The corrected phase image is input to a 
complex image producing section 720. The complex 
image producing section 720 is also supplied with the 
power image from the power image producing section 
702. The complex image producing section 720 pro- 
duces a complex image based on the phase image and 
the power image. The real and imaginary parts of the 
complex image can be obtained as the cosine and 
since, respectively, of the power image data. It should 
be noted that the angle for use in the cosine and sine 
calculation is a phase angle. 

[0064] The complex image is input to a phase mod- 
ification section 722, which is an embodiment of the 
phase-modification means of the present invention. The 
phase modification section 722 modifies the phase of 
the complex image based on a histogram input from a 
histogram generating section 724. 
[0065] The histogram generating section 724 gen- 
erates a histogram of the phase of the complex image. 
Thus, a histogram as exemplarily shown in Fig. 10 is 
obtained. As shown, the histogram has two peaks. A 
portion having its peak at a phase of zero corresponds 
to water, and a portion having its peak near a phase of 
jc/2 corresponds to fat. 

[0066] Generally, since fat contains water in the 
same voxel, its peak in the phase histogram shifts from 
tc/2. The reason is that a vector of fat signal F is the sum 
of a fat vector f and a water vector w, as exemplarily 
shown in Fig. 1 1 . The phase modification section 722 
modifies the phase of the fat image according to a peak 
shift 5 of the histogram of fat. This provides the phase 
difference 7i/2 between the water and fat images. 
[0067] The phase-modified complex image is input 
to a water/fat separation section 726, which is an 
embodiment of the image-producing means of the 
present invention. The water/fat separation section 726 
produces a water image using the real part of the com- 



plex image and produces a fat image using the imagi- 
nary part. The water image produced is stored in a 
water image memory 728, and the fat image in a fat 
image memory 730. 

5 [0068] Thus, according to the present invention, 
information representing static magnetic field inhomo- 
geneity is acquired from an image obtained by one-time 
imaging operation, a phase error due to the static mag- 
netic field inhomogeneity is corrected, and a water 

10 image and a fat image are separated using the image 
after being corrected for the phase error, therefore ena- 
bling imaging with water and fat separated to be per- 
formed very efficiently. 

is Claims 

1. A magnetic field inhomogeneity measurement 
method comprising the steps of: 
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obtaining a tomographic image of a subject 
within a static magnetic field space using mag- 
netic resonance, said tomographic image hav- 
ing a phase difference of 27i/n (n > 2) between 
pixel data of water and fat; 

multiplying the phase of pixel data of said tom- 
ographic image by n; 

correcting wraparound caused by said phase 
multiplied by n exceeding a range between ±jc; 
and 

multiplying said wraparound-corrected phase 
by 1/n to obtain a phase representing magnetic 
field inhomogeneity in said static magnetic field 
space. 

2. A magnetic field inhomogeneity measurement 
apparatus comprising: 

imaging means for obtaining a tomographic 
image of a subject within a static magnetic field 
space using magnetic resonance, said tomo- 
graphic image having a phase difference of 
2jc/n (n > 2) between pixel data of water and fat; 

phase-multiplier means for multiplying the 
phase of pixel data of said tomographic image 
by n; 

wraparound-correction means for correcting 
wraparound caused by said phase multiplied 
by n exceeding a range between ±n\ and 

magnetic field inhomogeneity-measurement 
means for multiplying said wraparound-cor- 
rected phase by 1/n to obtain a phase repre- 
senting magnetic field inhomogeneity in said 
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static magnetic field space. 

3. A phase correction method comprising the steps of: 

obtaining a tomographic image of a subject s 
within a static magnetic field space using mag- 
netic resonance, said tomographic image hav- 
ing a phase difference of 2n/n (n > 2) between 
pixel data of water and fat; 

w 

multiplying the phase of pixel data of said tom- 
ographic image by n; 

correcting wraparound caused by said phase 
multiplied by n exceeding a range between ±n\ 15 

multiplying said wraparound-corrected phase 
by 1/n to obtain a phase representing magnetic 
field inhomogeneity in said static magnetic field 
space; and 20 

correcting the phase of pixel data of said tomo- 
graphic image using said phase obtained. 

4. A phase correction apparatus comprising: 25 

imaging means for obtaining a tomographic 
image of a subject within a static magnetic field 
space using magnetic resonance, said tomo- 
graphic image having a phase difference of 30 
2n/n (n > 2) between pixel data of water and 
fat; 

phase-multiplier means for multiplying the 
phase of pixel data of said tomographic image 35 
by n; 

wraparound-correction means for correcting 
wraparound caused by said phase multiplied 
by n exceeding a range between ±rc; 40 

magnetic field in homogeneity-measurement 
means for multiplying said wraparound-cor- 
rected phase by 1/n to obtain a phase repre- 
senting magnetic field inhomogeneity in said 45 
static magnetic field space; and 

phase-correction means for correcting the 
phase of pixel data of said tomographic image 
using said phase obtained. 50 

5. A magnetic resonance imaging apparatus compris- 
ing: 

imaging means for obtaining a tomographic 55 
image of a subject within a static magnetic field 
space using magnetic resonance, said tomo- 
graphic image having a phase difference of 



2n/n (n > 2) between pixel data of water and fat; 

phase-multiplier means for multiplying the 
phase of pixel data of said tomographic image 
by n; 

wraparound-correction means for correcting 
wraparound caused by said phase multiplied 
by n exceeding a range between dbc; 

magnetic field inhomogeneity-measurement 
means for multiplying said wraparound-cor- 
rected phase by 1/n to obtain a phase repre- 
senting magnetic field inhomogeneity in said 
static magnetic field space; 

phase-correction means for correcting the 
phase of pixel data of said tomographic image 
using said phase obtained; and 

image-producing means for separately produc- 
ing a water image and a fat image using the 
phase difference in said phase-corrected pixel 
data. 

6. A magnetic resonance imaging apparatus as 
defined in claim 9, further comprising: 

phase-modification means for modifying the 
phase difference in said phase-corrected pixel 
data based on a histogram of the phase of said 
pixel data, 

wherein said image-producing means sepa- 
rately produces a water image and a fat image 
using said modified phase difference. 

7. A magnetic resonance imaging method comprising 
the steps of: 

obtaining a tomographic image of a subject 
within a static magnetic field space using mag- 
netic resonance, said tomographic image hav- 
ing a phase difference of 27c/n (n > 2) between 
pixel data of water and fat; 

multiplying the phase of pixel data of said tom- 
ographic image by n; 

correcting wraparound caused by said phase 
multiplied by n exceeding a range between :bc; 

multiplying said wraparound-corrected phase 
by 1/n to obtain a phase representing magnetic 
field inhomogeneity in said static magnetic field 
space; 

correcting the phase of pixel data of said tomo- 
graphic image using said phase obtained; and 
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separately producing a water image and a fat 
image using the phase difference in said 
phase-corrected pixel data. 



8. A magnetic resonance imaging method as defined 5 
in claim 7, further comprising the steps of: modify- 
ing the phase difference in said phase-corrected 
pixel data based on a histogram of the phase of 
said pixel data; and separately producing a water 
image and a fat image using said modified phase w 
difference. 

9. A method or apparatus as defined in any preceding 
claim, wherein n is equal to 4. 
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